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Abstract 

The difference between the number of lattice points in a disk of 
radius y/i/2n and the area of the disk t/4ir is equal to the error in 
the Weyl asymptotic estimate for the eigenvalue counting function of 
the Laplacian on the standard flat torus. We give a sharp asymptotic 
expression for the average value of the difference over the interval 
< t < R. We obtain similar results for families of ellipses. We 
also obtain relations to the eigenvalue counting function for the Klein 
bottle and projective plane. 

DD 

1 The simplest case 

Consider the standard flat torus [0, 1] x [0, 1] with boundaries identified. 
The eigenfunctions of the Laplacian are e 2mn ' x for n G Z 2 with eigenvalues 
(27r) 2 |n| 2 , so the eigenvalue counting function is 

N(t) = 4f (n G Z d : \n\ < v^/27r| , (1.1) 

the number of lattice points inside the disk B^q^ 2n of radius y/i/2ii about the 
origin. To first approximation N(t) is the area of the disk t/4ir, and this is 
exactly the Weyl asymptotic law. The problem of estimating the difference 

D(t) = N(t) - ± (1.2) 

is notoriously difficult (conjectured to be 0(t 1/,4+e ) for every e > 0). Here 
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Figure 1: D(t) Figure 2: t l ^D{t) 
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Figure 3: A(t) Figure 4: t l/A A{t) 

we study the simpler problem of approximating the average value 

i r R 

A(R) = - / D(t) dt. (1.3) 

Note that we are not taking the absolute value of D(t) in the average, so 
we may exploit the cancellation from regions where N(t) is greater than and 
less than t/An. We will show that A(R) = 0(R~ 1 ^ 4 ) as R — > oo, and more 
precisely 

A(R) = g(R~ l/2 )R~ l/4 + 0(R~ 3/4 ) as R -)• oo (1.4) 

where g(-R) is an explicit uniformly almost periodic function of mean value 
zero. Somewhat different but related ideas are given in Bleher j2j [3] . The 
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following Lemma is well-known (see jl], p. 74), but we include the proof for 
the convenience of the reader. 



Lemma 1. We have 



;i.5) 



where n = (ni,n 2 ) is a variable in Z 2 , and J 2 denotes the Bessel function. 
The series in (1.5) converges uniformly and absolutely. 

Proof. Let Xt denote the characteristic function of the ball B^.^ 2n . It is 
well-known that 

z = 



Xt{z) 



2tt\z\ ' 

4-7T 



Following standard methods (see [7] or [1]) we apply the Poisson summation 
formula to 

1 



R 



Xt * ^5 dt, 



(1.7) 



where ips is a smooth approximate identity. The ip$ convolution makes Frj 
smooth, but eventually we will let 5 — > 0. Note that 



The Poisson summation formula gives 

neZ 2 n&? 



;i.8) 



;i.9) 



V- 



R Jo 

R 



Xt(n)if)(5n) dt 



by Ql . 6p . Combining (|1.8|) and (jl.9 ) yields 

'* Vt 
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A(R) = \im\2l f 



Ji(\n\Vt)dtii>(6ri). 



[1.10) 



Now we use the property of Bessel functions Q6J) 

/ s a+1 J a (s)ds = R a+1 J a+l {R) (1.11) 
Jo 

for a = 1, together with the change of variables s = \n\y/i, to evaluate the 



integral in (1.10) 



R l n 2u\n 



Ji(\n\Vi)dt 



R J 7r|n| 4 

-^U\n\R), 
7r \n r 



J i(s) ds 



(1.12) 



and substitute this into (1.10) to obtain 



A(R) = lim ——J 2 (\n\\^RU(5n). 



;i.i3) 



The estimate J2(\n\yR) = Q( .11/^1/4 ) shows the convergence of the sum in 



(1.13) without the term ip(5n), so we can take the limit in (1.13) and obtain 

(TW. □ 



Theorem 2. Consider the uniformly almost periodic function with mean 
value zero 



9K X ) 



71 



3/2 



n| 5//2 cos 1 //!./• — 
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= g{R 1/2 )R~ 1/4 + 0(iT 3/4 ) as R -> 00. 



;i.i4) 



;i.i5) 



More generally, there exists a sequence of uniformly almost periodic functions 
gi, 02, ... witt gi = g such that for any n, 



A{R) = Y,9AR 1/2 )R l '"+0{R 
i=i 



1 n , 

' 4 2 



(1.16) 



Proof. We use the well-known asymptotic expression for Bessel functions 
"2 



JJx) 



1/2 cos ^x - \xtx - ^\ + 0(x 3/2 ) as x ->■ 00. (1.17) 
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Figure 5: g{\fi) 



Figure 6: t l ^A{t) - g(s/i) 



When a = 2 this is 

Ux) 



x 1//2 COS I x — 



x-l)+0(x~V% 



and we substitute this into (1.5) with x — \n\vR to obtain 



(1.18) 



(1.19) 



It is easy to see that the remainder term in (1.19) is 0(i? -3 / 4 ), so (1.19) yields 
(1.15). To obtain the more refined asymptotic expression ( 1.16[ ) we use the 
known more refined asymptotic expansion for Bessel functions (see [B]). In 
particular we note that it is possible to obtain explicit series expansions of 
the functions g^ for example, 



g-A x ) 



15\/2 

n^0 



n 



-7/2, 



7T 



sin \n\x — 



;i.2o) 



□ 



It is also reasonable to consider the function iV((27rr) 2 ) that counts the 
number of lattice points inside the ball B r of radius r, the difference D((27rr) 2 ) 
N((2nr) 2 ) — nr 2 , and the average with respect to the radius variable 



A(R) = - / D((2nr) 2 ) dr. 

R Jo 



(1.21) 
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This is a different average, but a change of variable shows that 

MR) = ^ R ( ttzW / ^ (*)^ 




;i.22) 



Since most of the contribution to the integral occurs for values of t near 
2ttR 2 , we see that A(R) has the same asymptotics as ^^A(2tiR 2 ). 

In Figure [l] we show the graph of D(t) and in Figure [2] the graph of 
t~ 1 / /L D(t). This illustrates the rough t 1 / 4 growth rate of D(t). In Figure 3 
we show the graph of A(t), and Figure [4] the graph of t^Aft). Figure 5 
shows the graph of g(Vi), which is almost identical to Figure kl for large t. 
Figure 6] shows the difference of t^^A^) and g(y/i), and Figure [7] shows this 
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difference multiplied by t x l 2 . Figure 8] shows the graph of A(t). Figure [9] 
shows the graph of -^^A(2nt 2 ), which agrees with Figure [8| for large £, and 



Figure 10 shows the difference. For more data see the website [3] 



2 The general case 

Consider the general flat 2-dimensional torus, M 2 /£ for some lattice C. The 
eigenfunctions of the Laplacian (restriction of the standard IR 2 Laplacian) 
have the form e 2mx for £ in the dual lattice with eigenvalues (27r) 2 £ • £. 
By diagonalizing the quadratic form £ • £ on C we can find an orthonormal 
basis t>i,t>2 in IR 2 and positive constants a\, a 2 , such that the eigenvalues are 

<-> 2 ((^) 2 + (^) 2 ) — *■ 

Thus the eigenvalue counting function is 

((^H^)T^ t/ar >- (21) 

In place of disks we consider the family of ellipses 

Of course iV(t) is just the number of lattice points in E t , and the volume of 
E t is 2122*. Again we write £>(t) = N(t) - for the difference and define 



the average A(i?) by (1.3). The analog of (1.6) is 



- : Ji(^J(aiZ ■ vi) 2 + (a 2 z ■ v 2 ) 2 Vi) z^O 



\a^t Z=Q 

(2.3) 

Lemma 3. We have 

A(R) = —r ^-Tl WiMV ( a l n - U l) 2 + ( a 2 n - V2) 2 VR), 

^n{am ■ Vl ) 2 + [a 2 n ■ v 2 ) 2 ] 

(2-4) 

the series converging uniformly and absolutely. 
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sooo inooo 



Figure 11: g(\ / t),a 1 = 2,a 2 = 1/2 



Figure 12: t l l A A(t) - g(Vi) 



Proof. The same proof as for Lemma [TJ with (2.3) used in place of (1.6). □ 



Theorem 4. The asymptotic expansions (1.15) and (1.16) hold, where now 

V2 



9{x) 



3/2 



7T 



■ COS 



aia 2 2^ (( a i n ' ^i) 2 + ( a 2"- • v 2 ) 2 ) 



2\"5/4 



\( ai n ■ v,) 2 + (a 2 n ■ v 2 ) 2 ) 1/2 x - ^) . (2.5) 
Proof. Same as for Theorem |2j using Lemma [3] in place of Lemma [T} □ 



See Figure 
difference t 
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for g(y/i) with a\ = 2 and a 2 = 1/2 and Figure 12 for the 
^A(t) - g(y/i) for the same family of ellipses. Plots for different 
families of ellipses are available on the website [5]. 



3 The Klein bottle and projective plane 

If we identify the vertical boundaries of the square directly, and the horizontal 
boundaries with reflection, we obtain the standard flat Klein bottle KB. In 
terms of functions defined on the square, we are imposing the boundary 
conditions u(0,y) = u(l,y) and u(x, 0) = u(l — x,l) in order to have a 
function on KB. We may cover KB by the rectangular torus [0,1] x [0,2] 
with the identities 

u(x + l,y) = u(x,y) 
u(l -x,y + l) = u(x,y) 



(3.1) 
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describing the lifts of functions on KB to M 2 . The eigenfunctions of the 
Laplacian on KB lift to eigenfunctions on the rectangular torus, and so 
are linear combinations of functions of the form e 2m x+ 2y) with eigenvalue 
(2vr) 2 (j 2 + (§) 2 ). Now we observe that e *rf(i(i-*)+f (v+i)) = (_i) 
Thus there are two families of eigenfunctions 

e 2m ? y for k even (corresponding to j = 0), and (3.2) 
e 27ri (^+f*0 + (_ 1 jfc e 2m(-ix+| w ) for ^ > Q (3 3) 



We can therefore see that the eigenvalue function Nkb is close to one half 
the counting function iV Tl 2 for the [0, 1] x [0,2] torus. 

Theorem 5. N KB (t) = |iV Tl2 (t) ± \. 

Proof. iV" Tl 2 (t) counts all integers j,k such that j 2 + (|) 2 < t^w- When 
j ^ the pair ±j contributes just a single eigenvalue to A^kb(*)- When 
j = we count all k such that < — in iV Tlj2 (t), but just the even values 
of k in N KB (t), and even : \k\ < ^} = \{k : |A;| < ^} ± |. □ 

It is interesting to compare the Klein bottle with the projective plane 
(PP) obtained from [0, 1] x [0, 1] by identifying both sets of boundary edges 
with reflections. Functions on PP lift to 1R 2 with the identities 

w(l - x, y + 1) = u{x, y) 
u(x + 1,1 - y) = u(x,y) 

and the torus [0,2] x [0,2] is a four-fold covering of PP. However, while it is 
possible to pull back the standard Laplacian to PP, the pairs {(0, 0), (1, 1)} 
and {(0, 1), (1, 0)} of identified points on PP are singularities (cone points 
with total angle 7r) with respect to the otherwise flat metric. 

Reasoning as in the KB example, we know that eigenfunctions of the 
Laplacian on PP must be linear combinations of the functions e 2m ^ x+ 2 y ^ 



with eigenvalue (27r) 2 ((|) 2 + (o) 2 )- Imposing the conditions (3.4) leads to 
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four families of eigenfunctions: 

constants (corresponding to j = and k = 0) (3.5) 
e 2l 4y + e - 2m l^ f or > o even (corresponding to j = but fc ^ 0) (3.6) 
e 27ri|x _j_ e -2-Ki\x £ Qr j > o even (corresponding to = but j ^ 0) (3.7) 

for j > and jfc > 0. (3.8) 

This leads to the identity 

iVpp(t) = \N T2a (t) + \±\ = \N Thl (4t) + i ± \. (3.9) 

Similar results hold for KB and PP constructed from the tori considered 
in section [2] Related questions in the context of fractal Laplacians with the 
Sierpinski carpet replacing the square are discussed in [T]. 
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